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1. INTRODUCTION 

In this study we present a series of Large Eddy 
Simulations (LES) employing the Super-Droplet 
Method (SDM) for represe nting aerosol , cloud 



and warm-ra in microphysics flShima 2008 Shima 
|et al.| |2009[ ). SDM is a particle-based and prob 
abilistic Monte-Carlo type model. The particle- 
based formulation helps to overcome the prob- 
lem of parameterisation of processes occurring 
at single-particle scale (micro- to millimetres) as 
source/sink terms in the LES equations solved on 
grids with cell dimensions of the order of tens of 
metres. The aim of this paper is to showcase the 
capabilities and point out the limitations of SDM 
in context of simulation of a field of precipitating 
clouds. 

2. THE SUPER-DROPLET METHOD 

The framework of the method consists of 
two mutually coupled simulation components: 
(i) a fluid flow solver computing (in an Eule- 
rian sense) the evolution of fluid velocity field 
and evolution of the thermodynamic scalar quan- 
tities, and (ii) a particle-tracking logic comput- 
ing (in a Lagrangian sense) evolution of physical 
coordinates and physicochemical properties of a 
population of particles. The coupling between 
the Eulerian and Lagrangian components is bi- 
directional. The Lagrangian component feeds on 
the fluid velocity field in order to update the po- 
sitions of super-droplets, and the thermodynamic 
fields to compute the condensational growth or 
evaporation rates. The Eulerian component feeds 
on the water vapour and heat source/sink rates 
resulting from condensation and evaporation of 
water on the particles. Although mutually cou- 
pled, the Eulerian and Lagrangian computations 
are decomposed, in the sense that they happen 
sequentially and the state vectors of one com- 
ponent are constant from the standpoint of the 
other (within a time-step). 

The model does not differentiate between 
aerosol particles, cloud droplets, drizzle or rain 
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drops. Each particle in the model (referred to 
as super-droplet) represents a multiplicity of real- 
world particles of the same size and of the same 
chemical composition. From the standpoint of 
the Lagrangian component, the super-droplets 
are subject to three processes besides advec- 
tion: (i) gravitational settling, (ii) condensational 
growth/evaporation and (iii) col lisiona I growth. 
Consequently, the model covers representation of 
such cloud-microphysical processes as: cloud con- 
densation nuclei (CCN) activation; drizzle forma- 
tion by collisions of cloud droplets (autoconver- 
sion); accretion of cloud droplets by drizzle drops 
and raindrops, as well as coalescence of these 
larger hydrometeors (self-collection); and precip- 
itation of drizzle and rain including aerosol wet 
deposition. Analogous particle-based techniques 
were recently used in context of si mulations of 
precip itation-formin g clouds e.g. by |Andrejczuk| 
|et a I . | Q2010D and [Franke and Raaschl ( |2010| ); 
what distinguishes SUM however, is the prob- 
abilistic Monte-Carlo type representation of the 
particle coalescence process. Probabilistic repre- 
sentation of particle coalescence has also been 
recently used in atmospheric modelling studies; 
however, to authors' knowledge, none of these 
models w ere bi-directiona ll y cou pled with LES - 
see e.g. |Jensen and Lee (2008) for description 
of an adiabati c parcel mo del with Monte-Carlo 
coalescence or |Riemer et al. ( |2009[ ) for descrip- 
tion of particle-resolved aerosol transport model 
with Monte-Carlo coalescence but working off- 
line from the flow dynamics). 

3. MODELLING SET-UP 

All simulations in the present study were car- 
ried out using the Nagoya Universit y Cloud- 
Resol ving Storm Simulator (CReSS: |Tsuboki| 
2008, chapter 9.2 and references therein). CReSS 
is a non-hydrostatic compressible flow LES-type 
solver. Two types of moist processes represen- 
tations were used: the SDM and a bulk mi- 
crophysics model (Kess ler-type parameterisation 



implemented following |Klemp and Wilhelmson 
1978 section 2b). The simulations are carried 
out using a set-up b ased on the RICO com- 
posite case defined in van Zanten et al. (2011) 
and corresponding to atmospheric state measured 
and modelled in context of the Rain in Cumulus 



over Ocean (RICO) field project flRauber et al. 
2007). The set-up defines initial profiles of po 



Table 1: List of model runs discussed in the paper. The run label denotes whether bulk (blk) or SDM (sdm) microphysics 
was used, as well as which grid resolution (coarse, middle or high) and super-droplet number density was chosen. Coarse 
resolution corresponds to a quarter of the domain from the original RICO set-up (i.e. grid box size of 100x100x40 m 
with 64x64x100 grid points); the middle and high resolutions denote settings resulting in halved and quartered grid 
box dimensions, respectively (with the domain size kept constant). For each simulation there are five time-steps defined: 
long and short time-step of the Eulerian component (the short one used for sound-wave terms), the time-step used for 
integrating the condensational growth/evaporation equation, the time-step used for solving collisional growth using the 
Monte-Carlo scheme, and the time-step for integration of particle motion equations. 



run label 


grid 


dx=dy 


dz 


time-steps [s] 


sd density 


cm" 3 ] 


blk-coarse 


64 x 64 x 100 


100m 


40m 


1.00/0.100 n/a 


n/a 




sdm-coarse-8 


64 x 64 x 100 


100m 


40m 


1.00/0.100/0.25/1.0/1.0 


2.0x10 


-li 


sdm-coarse-32 


64 x 64 x 100 


100m 


40m 


1.00/0.100/0.25/1.0/1.0 


8.0x10 


-li 


sdm-coarse-128 


64 x 64 x 100 


100m 


40m 


1.00/0.100/0.25/1.0/1.0 


3.2x10 


-10 


sdm-coarse-512 


64 x 64 x 100 


100m 


40m 


1.00/0.100/0.25/1.0/1.0 


1.3x10 


-09 


sdm-middle-8 


128x128x200 


50m 


20m 


0.50/0.050/0.25/1.0/1.0 


1.6x10 


-10 


sdm-middle-32 


128x128x200 


50m 


20m 


0.50/0.050/0.25/1.0/1.0 


6.4x10 


-10 


sdm-middle-128 


128x128x200 


50m 


20m 


0.50/0.050/0.25/1.0/1.0 


2.6x10 


-09 


sdm-high-8 


256x256x400 


25m 


10m 


0.25/0.025/0.25/1.0/0.5 


1.3x10 


-09 



tential temperature, specific humidity and wind 
characteristics for the trade-wind boundary layer 
follow ing RICO observations (see |Nuijens et aL~| 
[20091 ). Th 's set -up was previously used in sev- 



eral m o delling studies including: ([Stevens anc 
Seifert, 2008; Jiang eTatj |2 009l Ibeitert et a l 
2010; Matheou et al.TMT 



Grabowski et a I 



201ip . Here the only exception from the origina 
set-up is the domain size - a quarter of the origi- 
nal domain was used (i.e. 6.4x6.4x4 km instead 
of 12.8x12.8x4 km), and the grid cell sizes - 
several settings were used. Tests with full domain 
size (not discussed herein) revealed that quarter- 
ing the domain enlarges the fluctuations in time 
of cloud macroscopic properties; however, not to 
a level significant for the presented discussion. 

For simulations using SDM, the initial coordi- 
nates and sizes of particles are chosen randomly 
with uniform distribution in physical space and 
bi-modal lognormal distribution in particle-size 
space (size spectrum as specified in van Zanten 



et al 



2011 section 2.2.3). All particles are ini 
tially put in equilibrium with ambient humidity 
assuming that all are composed of ammonium 
sulphate solution. The mean number of super- 
droplets per LES grid box varied from 8 to 512 in 
different simulations (see Table Tj). | n all calcu- 
lations discussed in the paper trie probability of 
collis i ons a nd coalescence was defined following 



Hall (1980 



hence no e 



section 3d, and references therein) 
fects of small-scale turbulence on the 



drop collision efficiency was taken into account. 

4. SIMULATION RESULTS 

A list of model runs and their correspondin 
labels used in the paper is given in Table 
The SDM simulations were carried with differ 
ent grid resolutions and different super-droplet 
number densities. The table also lists different 
time-steps used in the simulations - in SDM the 
time-steps used for the Eulerian and Lagrangian 
computations differ as the corresponding numer- 
ical stability constraints differ. 



Fol lowing the methodology of van Zanten et ah] 
( |2011[ ) the analyses are restricted to the last tour 



hours of the 24-hour model runs. Presented re- 
sults are based on the LES grid data output ev- 
ery minute simulated time and super-droplet data 
(particle positions and sizes) output every 10 min- 
utes simulated time. 

4.1 CLOUD MACROSTRUCTURE 

Figure [Tj presents time-series of scalar quanti- 
ties characterising domain-wide cloud field fea- 
tures: cloud cover (cc), liquid water path (LWP), 
rain water path (RVVP) and cloud top height 
(zct). All quantities are labelled a nd defined fol- 
lowing |van Zanten et aT. (2011 section 2.4.2 
and Table 4 therein, ana caption of Fij 
herein). Presented time-series cover the first a 
the last four hours of simulations corresponding 
to the spin-up, and the relatively steady-state 
stages of the simulation, respectively. The figure 
is intended for comparison with analogous plots 



up: Fig. 2 in 


Stevens and Seifert (2008), Fig. 5 in 


Seifert et al. 


(2010), Fig. 


3 in van Zanten et al. 


(2011), Figs 1-3,7 in 


Matheou et al.| Q2011 ) and 


Fig. 6 in 


Grabowski et al. 


(20111). 



The cloud-cover (cc) plots reveal considerable 
dependence of cc on the choice of grid resolu- 
tion, with a significant increase of cloud cover 
for higher-resolution simulations. This is, at least 
partially, caused by the definition of cc which in- 
cludes an arbitrary threshold value for the cloud 
water mixing ratio (e.g. in a hypothetic 2x2 
chessboard water content distribution with cc of 
50%, downsampling to a single grid point could 
result in a zero cc). Furthermore, since the forma- 
tion of convective clouds is triggered by vertical 
air motion and since refinement of the grid resolu- 
tion helps to resolve in more detail the dynamics, 
an increase in the number of (smalle r) clouds may 
be expected ( |Matheou et al.||201l| Sect. 3.2). 

Liquid water path (LWP) plots show that the 
increase of cc with increased resolution is corre- 
lated with increase of LWP. The values of LWP 
obtained in middle-resolution SDM simulations 
(blue lines) fall within the range of v alues ob 
tained with other LES from the study of van Zan 




Rgm^^Tim^^ene^o^joudmacro^ 



isTeTinecf^is^n'e^r3cTlo7T^f^oTLm*Tns^<vrH grid cell (a grid cell is defined 

as cloudy if the cloud water mixing ratio exceeds 0.01 g kg" 1 ); (ii) liquid water path " LWP" is the mean over all columns of 
column-integrated liquid water content (density of both cloud and rain water); (iii) rain water path " RWP" is calculated in the 
same manner using rain water only; (iv) the cloud top height "zct" is the height of the top of the highest cloudy grid cell. For 
SDivi simulations the cioud and rain water mixing ratios are diagnosed by summing over the popuiation of particies with radius 
smaller and larger than 4U um. resoectiveiv. in the foreground there are clots deciding data from the nine simulations listed in 



therein ): ^ cloud cover "cc" 



I a ble fl| In the background (plotted with thin grey lines) there are the resu lts from Lbb simulations described in van Zanten et at. 
^2o"^fjaa*ta" obtained at: http://knmi.nl/samenw/rico/ ^^ee^ection |4. 1| for discussion - 
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Figure 2: Height-resolved statistics of cloud droplet spectrum param eters derived from SDM simulation data. All plots in a row 



share the same y-axis representing height above sea level. See section pT^|tor discussion and plot construction method description 



□ 



ten et a I 



( 2011| grey lines in Fig. JI herein) for 
most ot the simulated time. The high- resolution 
run (green line) features highest fluctuations of 
values of both LWP and RWP which likely cor- 
responds to the low super-droplet density, and 
hence the low spectral resolution of the micro- 
physics model. 

The rain water path was one of the quanti- 
ties that varied most from one simulation to an- 



other in the study of van Zanten et a\] Q2011| 
Sect. 3.2). Similarly, the hereby presented re- 
sults exhibit both significant variations in RWP 
with simulation time as well as significant depar- 
tures from one simulation to another (i.e. when 
changing grid resolution or super-droplet density). 
A notable feature is the high correlation of RWP 
with both LWP and cc. 

The plots of cloud-top evolution with time 
clearly show that the SDM-simulated cloud field 
is shallower than the one obtained with bulk mi- 
crophysics or in other LES. While the employed 
definition of cloud-top height does not allow di- 
rect comparisons with measurements, arguably 
the SDM-predicted heights for middle- and high- 
resolution seem more com parable to those ob 



served during RICO ( e.g. Rauber et al. 2007 



significant inhibition of convection in all coarse 
resolution runs suggest that the 100 x 100 me- 
ter resolution is not enough when coupled with 
a Lagrangian model which, among other factors, 
relies on LES-predicted supersaturation field di- 
rectly coupled with the vertical motions. 

4.2 CLOUD MICROSTRUCTURE 



Genkova et al. 2007| Fig. 1). In general, the 



Figures 2|3 and [4] present height-resolved 
statistics ot the vertical velocity w, the super- 
saturation S, cloud droplet concentration CDNC, 
droplet effective radius r e ff=<r 3 > / <r 2 >, liq- 
uid water content LWC, the cubed ratio of mean 
volume radius to effective radius k=<r 3 > / r 3 ^., 

and the standard deviation of cloud droplet ra- 
dius cry. The plots are intended for compari- 
son w ith the analysis presented in Arabas et al. 
(2009, Figs. 1 and 2) where the data from aircraft 
measurements during the RICO campaign using 
the Fast-FSSP optic a l clou d droplet spectrome- 
ter (Brenguie r et al.| |1998[ ) were analysed. The 
herein analysis ot bUM simulation data is con- 
strained to in-cloud regions defined as the grid 
boxes having CDNC> 20 cm" 3 where CDNC is 
derived by summing over the super-droplets rep- 
resenting particles of radius between 1 and 24 mi- 




crometers. The choices of the CDNC threshold 
and the spectral range correspond to those char- 
acteristic of the Fast-FSSP probe. 

Plot construction method was chosen follow- 



ing the methodology of Arabas et al. (2009). For 
each level of the model grid and each plotted pa- 
rameter a list of values matching the in-cloud cri- 
terion is constructed, sorted and linearly interpo- 
lated to find the 5 th , 25 th , 45 th , 55 th , 75 th and 
95 th percentiles. The lists are constructed from 
the LES-grid values (w, S) or super-droplet statis- 
tics calculated for each grid cell (CDNC, r e ff, 
LWC, k and oy). The sample volumes are there- 
fore defined by grid cell sizes and are of the order 
of 10 4 - 10 5 rn 3 , while the sample volume for the 
Hz Fast-FSSP dataset used in lArabas et al.l 



(2009) is of the order of 10 



-6 



771. 



a thing one 

has to bare in mind. The vertical extent of the 
measurement sample volume of ca. 10 m is com- 
parable to the grid cell size, though. Levels with 
less than 0.1 • N max data-points, where N max is 
the number of data-points at the level with the 
highest number of data-points, are discarded from 



the analysis. Finally, the 5 - 95 percentile, 
the interquartile, and the 45 th - 55 th percentile 
ranges are plotted as a function of height using 
red, green and blue bars, respectively. The pro- 



files composed of the 45 th 55 th percentile range 
bars are referred to as median profiles hereinafter. 

VERTICAL VELOCITY Profiles of the verti 
cal velocity generally show a gradual increase of 
the spread of values with the increasing grid reso- 
lution. It is consistent with the fact that the grid 
acts as a spatial filter in the LES. The median 
profiles resemble each other regardless of the res- 
olution suggesting that the macroscale dynamics 
of the simulated cloud layer are robust to both 
the LES grid choice and the super-droplet den- 
sity choice (however, the differences in cloud-top 
heights plotted in Fig. [l] and evident in Figs.[2]-[4] 
reveal that this robustness is limited). The slight 
increase of the vertical velocities with height visi- 
ble in some of the profil es, and noted i n the RICO 



observations analysis (Gerber et al. 



2008 Ta 



ble. 2), may correspond to the tact that only the 
more vigorous updraft regions were able to pro- 
duce clouds reaching the upper part of the cloud 
field. 

SUPERSATU RATION The profiles of super- 
saturation, especially for the high and middle res- 
olution, show the characteristic cloud-base max- 
ima. This confirms that the model formulation, 
and the time resolution used allows to capture 
the influence of CCN activation kinetics on the 
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Figure 4: Same as Fig.[2]for model runs: sdm-middle-8, sdm-middle-32 and sdm-middle-128. 



evolution of supersaturation. Both the values of 
maximal supersaturation (ca. 1%), as well as a 
roughly estimated supersaturation relaxation time 
(ca. 100s assuming «lm/s vertical velocity, and 
an m 100m height scale over which the supersat- 
uration falls off to an asymptotic value) corre- 
spond to those reported in modelling studies that 
employed supersaturation- predict i ng models (e.g. 
M o r r i so n a nd G ra b ows k i | 2008| |Khvorostyanov| 



and Curry "2008). Comparison of the supersatu- 
ration prediction in the model with measurements 
is not viable as direct measurements of supersat- 
uration in clouds are v irtually unavailable ( |Ko- 



rolev and Mazin 2003). In none of the runs 
there is any other maximum of supersaturation 
visible along the profile. This suggests that the 
employed Lagrangian technique for representing 
water condensate inhibits the spurious production 
of cloud-edge supersatu ration inherent in Eulerian 



[20 
lik 



mode ls, (see Sec. 1 in Grabowski and Morrison 



2008, and references ther ein). The reason for itls 



ikely (as noted as well in |Andrejczuk et al~ 2008, 
par. 10) that such coupled Eulerian-Lagrangian 
approach does in fact cover representation of frac- 
tional cloudiness within a grid cell. 

DROPLET CONCENTRATION The range 
of drop concentration values obtained in the sim- 



ulations does roughly correspond to t he values ob- 
served during RIC O and presente d in|Arabas et al.| 
(2009 1 figure la), |Gerber et a I | Q2008j figure 4) 
and |Col6n-Robles et al.| ( |2006| figure 3). The dis- 
crepancies in the median values and the spread 
of CDNC among different model runs show, how- 
ever, that the prediction of drop concentration 
is sensitive to the super-droplet density (Fig. [4]). 
Furthermore, the higher concentrations obtained 
with low super-droplet densities (i.e. low spec- 
tral resolution) are closer to the observed values. 
The observed invariability of CDNC with height is 
robustly reproduced suggesting that the discrep- 
ancies are solely related to the treatment of CCN 
activation in the model. 

EFFECTIVE RADIUS The effective radius 
profiles are generally robust to both grid and 
super-droplet density choices, and they do re- 
semble the profiles observe d with the Fast-FSSP 
instrument during RICO QArabas et~aTj |2009 
fig. 2). The profiles show a gradual increase 
cloud droplet sizes from the cloud base up to the 
altitude of 700-800 metres where the median val- 
ues of r e ff reach 15-20 micrometres. Above, the 
profiles differ more from one model run to another 
but still the flattening of the r e fj profile slope is a 
robust feature. The reduced slope of the median 
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Figure 5: Comparison of model-predicted spectra for particles within size range of the 0AP-2DS instrument, within grid cells with 
ra in water mixing rati o q r > 0.001 g/k g and within an altitude range of 183 ± 100 m. Figure intended for comparison with Fig. 4 
in |Baker et al.| | [2009| |, see section |"4.3| herein for details. 



profiles of r e jf reflects (i) the chosen drop ra- 
dius range (1-24 fim - emulating the Fast-FSSP 
instrument range), (ii) the decreased efficiency, 
in terms of radius change, of the condensational 
growth for larger droplets, and (iii) the increased 
efficiency of drop collisions reached after the ini- 
tial conden sational growth sta ge (see e.g. sec- 
tion 3.1 in |Zhang et al. 2011, for discussion of 
influence of precipitation on effective radius pro- 
files). The spread of r e // values high above cloud 
base indicates presence of smaller droplets even 
near cloud top. The present analysis (and the 
corresponding choice of data output rate during 
simulations) does not provide an answer to the 
question of the origin of these small particles - 
e.g. if they were activated at cloud edges abov e 
cloud base or not (see e.g.|Slawinska et aL 2012[). 



LIQUID WATER CONTENT LWC profiles 
depicted in Figures |2]|4] show significant spread 
of values at a given level. This is in accord with 
RICO observations, and suggests mixing-induced 
dilution of cloud and the resultant decrease of 
water content (cf. the discu ssion of the adiabati c 
fraction profiles in Fig. 2 in Araba s et al.| |2009] ). 
In SDM the condensational growth and evapora- 
tion is computed using values of supersaturation 
interpolated to super-droplet positions, hence the 
mixing scenario is not homogeneous as different 
droplets are exposed to different supersaturations 
within a single grid cell. 

The decrease in LWC in the upper part of the 
cloud is correlated with an increase of rain water 
mixing ratio (not shown). In the lower part of 
the cloud, the median profiles of LWC obtained 



in the simulation with middle and high grid reso- 
lution show arguably reasonable agreement with 
LWC profi les derived from R I CO ob servation s de- 
picted in |van Zanten et a l. (2011 Fig. 8), |Ger- 
|ber et al.| ( |2008| Fig. 1) and |Abel and Shi pway 
(2007] Fig. 7) taking into account that the ref- 
erenced observations were taken with different 
instruments having different sampling rates and 
spectral ranges. 

PARAMETER K Most GCMs employ a pa- 
rameterisation of the effective radius when di- 
agnosing the cloud optical depth of clouds from 
predicted liquid water content and dro p number 
conce n tratio n (see e.g. section 2 in Brenguier 
et al. |2000[ and references therein). Typically, 
the cubed ratio of the mean volume radius to 
the effective radius fc=<r 3 > jr\^ (or a simi- 
larly defined scaling coefficient) is assumed to be 
constant or to depend solely on dro p concentra- 
tion (e.g. |Peng and LohmaTij |2003[ ). Value of 
k = 0.8 ± U.UY based on aircraft observations in 
maritime cl ouds using an FSSP-type instrument 
reported in Martin et al. ( |1994[ ) is often used in 
climate models! 

As shown in Figure [2] all simulations with the 
mean value of 8 super-droplet per LES grid cell 
predict values of k hardly different from unity. 
This means that the effective radius is in most 
cases equal to the mean volume radius, as it 
should be, taking into account that the 8 super- 
droplets (which may be thought of as size bins) 
are used to represent a wide spectral range of 
particles: aerosol, cloud, drizzle and rain; while 
the statistics presented in the plot cover a narrow 



spectral range. Consequently, cloud water is likely 
represented by a single super-droplet (bin) in each 
grid cell, a situation in which all of the cloud 
droplets would have the same size and k = 1 by 
definition. The values of the k parameter calcu- 
lated in the model runs with mean super-droplet 
density of 128 range from approximately 0.6 to 
0.9. There is no significant difference when the 
super-droplet density is increased to 512 suggest- 
ing that the density of 128 per cell is enough to 
resolve the relevant cloud droplet size spectrum 
features. 

DROP RADIUS STANDARD DEVIATION 

In consistency with the aforementioned behaviour 
of k for low super-droplet density, the drop ra- 
dius standard deviation a r hardly differs from zero 
when a mean of 8 super-droplets per cell is used, 
and it still goes down to zero at some levels with 
a mean of 32 super-droplets per grid cell. For all 
higher spectral resolution simulations the 5 th per- 
centile of a r is greater than zero at all levels, and 
the increase of mean super-droplet density from 
128 to 512 does not influence the profile shape. 
In general agreement with RICO observations, the 
standard deviation ranges from 1 to 6 microme- 
tres. The profiles of a r show a slight increase 
with height (best captured in the sdm-middle-128 
run); however, the observed inclination of the a r 
profiles resembles more the 95 th percentile profile 
derived from SDM simulations. The fact that the 
values of a r are larger than those obt ainable in 
adiabati c drop growth process (see e.g. |Yum and | 



Hudson| |2005| and references therein) suggests 
that the model does capture to some extent the 
mixing-induced broadening of the cloud droplet 
spectrum. However, since the mixing in SDM is 
limited to LES-resolved motions with character- 
istic length scales constrained by grid cell dimen- 
sions, the level of agreement with observations is, 
as expected, limited. 

4.3 PRECIPITATION MICROSTRUCTURE 

The analysis depicted in Figures [2l [3] and [4] 
concerned the cloud-droplet region oftne particle 
size spectrum. A comparison of model-predictec 
spectra for larger particles is depicted in Figure[5 
T he figure is i ntended for comparison with Fig. '■ 
in|Baker et al.|(|2009|) based on measurem ent data 



obtain ed" wit h the OAP-2DS instrument ( Lawson| 
et a 1. 1 12006"1) during RICO research flights. Dur- 
ing KlLU the OAP-2DS instrument was set to 
classify particles into 61 size bins spanning the 
2.5 M m - 1.5 mm size rang e in radius. In the anal- 
ysis of Baker et al .| (|2009 ) a mean size spectrum 
was derived from 26 1 spectra measured within 
rain-shafts below the cloud base at the altitude 
of about 183 metres (600 ft). In order to derive 
comparable quantities from the SDM simulation 
results, the super-droplets in each grid cell were 
classified into size bins of the same layout as used 
by the 0AP-2DS instrument, an altitude range of 
183 ± 100 m was chosen, and only grid cells with 
rain water mixing ratio q r > 0.001 g/kg were 



taken into account (q r being derived from sum- 
mation over super-droplets representing particles 
with radii greater than 40 ji tm). A comparison 
of Figure |5] with Fig. 4 from |Baker et al | ( [2009] ) 
(both plots share the same axis ranges) reveal 
that the model results, regardless of the grid- 
resolution or super-droplet density choice, show 
fair agreement with the measurement results for 
drop diameters greater than 0.1 mm. The spec- 
trum from the sdm-coarse-512 run (with the 
highest spectral resolution and the best Monte- 
Carlo sampling density) resembles most closely 
the measurements having the lowest concentra- 
tions in the 0.1-0.4 mm diameter range. All simu- 
lations disagree markedly with the measurements 
within the 10-20 /im diameter range. These mea- 
sured particles within t his range were " most likely 



deliquesced aerosols" ( |Baker et al. 2009) and 
since there are no aerosol sources in the model, 
this disagreement is still a plausible result. The 
drop breakup process was identified as another 
possible source of droplets smaller than 100 ^m, 
and this process is also not included in the present 
version of the model. Moreover, it was the lack 
of hydrometeors smaller than 100 /x in diameter 
that was considered as the primary high l ight o f 
the observations reported in Baker et al. (2009), 
and this feature of the spectrum is in tact hinted 
in the simulation results. 

5. SUMMARY 

The salient features that distinguish SDM from 
bulk and bin warm-rain microphysics models are: 

(i) diffusive error-free computational scheme for 
both condensational and collisional growth; 

(ii) particle spectrum representation allowing 
straightforward comparison with experimental 
data obtained with particle-counting instruments; 

(iii) persistence of arbitrary number of scalar 
quantities assigned to a super-droplet (e.g. chem- 
ical properties); 

(iv) scalability in terms of sampling error 
(i.e. super-droplet density); 

(v) parameterisation-free formulation of the key 
processes involved in cloud-aerosol interactions. 

The last feature, in particular, does not come 
without its challenges. Explicit treatment of 
aerosol microphysics implies taking care of their 
budget within the simulation domain, including 
modelling their sinks (wet and dry deposition) and 
sources, the latter not being represented in the 
present study. 

The other two processes not represented in the 
discussed simulations are: (i) the im pact of turbu- 
lence on drop collisions (see e.g. jDevenish et aTj 
|2012[ ), and (ii) the influence of drop breakup on 
the size spectrum of precipitation particles (see 
e.g. Villermaux and Bossal 2009). Yet, the ar- 



guably fair level ot agreement ot the simulation 
results reported herein and the previously pub- 
lished in-situ measurement results from RICO, 
suggest that the set of processes represented in 
the present set-up of SDM includes at least the 



key players involved in determining the features 
of the size spectra of cloud and precipitation par- 
ticles. One has to bare in mind that a direct 
comparison of RICO measurement data with LES 
results of the type presented herein is not possible 
due to different time- and space- scales associated 
with the model variables and the measurements, 
as well as due to the nature of the "compos- 
ite" modelling set-up. (typical atmospheric con- 
ditions, with the diurnal cycle neglected in partic- 
ular). 

The SDM modelling approach offers the unique 
possibility to mimic in the analysis the particle- 
counting process inherent in the principle of oper- 
ation of most airborne aerosol, cloud and precipi- 
tation probes. As a result, it becomes meaningful 
to analyse such model-predicted cloud droplet size 
spectrum parameters as e.g. k or a r , not taken 
into account in previous comparisons of RICO 



LES resu l ts with aircraft observat i ons ( Abel and 



Shipway| |2007| |van Zanten et al.| [20lip . What 
militates in favour ot pertinence ot the obtained 
results is that the present model employs fewer 
parameterisation (in comparison with bulk or bin 
models) and more basic principles to describe pro- 
cesses occurring at the microscale (e.g. descrip- 
tion of the Kohler curve shape as opposed to em- 
ployment of such parameters as the exponent in 
Twomey's formula for CCN activation parame- 
terisation in bin models, or the autoconversion 
threshold in Kessler-type bulk models). 

The key conclusion from the presented analysis 
is that while the convergence of the macroscopic 
cloud parameters do not seem to get any better 
than in the other LES simulations using the RICO 
set-up, the SDM is able to provide more detailed 
insight into cloud microstructure, and thus indi- 
rectly into its optical properties. 
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